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ABSTRACT 

We have obtained high dynamic range, good natural seeing i' images of BL Lacertae objects (BL 
Lacs) to search for the AGN host and thus constrain the source redshift. These objects are drawn 
from a sample of bright flat-spectrum radio sources that are either known (via recent Fermi LAT 
observations) gamma-ray emitters or similar sources that might be detected in continuing gamma- 
ray observations. All had spectroscopic confirmation as BL Lac sources, but no redshift solution. 
We detected hosts for 25/49 objects. As these galaxies have been argued to be standard candles, 
our measured host magnitudes provide redshift estimates (ranging from 0.2-1.0). Lower bounds are 
established on the redshifts of non-detections. The mean of the fit redshifts (and lower limits) is higher 
than those of spectroscopic solutions in the radio- and gamma-ray- loud parent samples, suggesting 
corrections may be needed for the luminosity function and evolution of these sources. 
Subject headings: galaxies : active - galaxies: BL Lacertae objects: general 
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1. INTRODUCTION 



BL Lac objects, being highly continuum-dominated, are a perennial problem for those wishing to study the evolution 
of AGN populations. The lack of visible broad lines and, even more, the low equivalent width of host absorption features 
makes redshift determinations extraordinarily difficult. Yet, since the conti nuum domination is a m anifestation of the 
£j ' good alignment of the relativistic jet outflow to the Earth line-of-sight, ([Urrv fc Padovariil 119951) knowledge of the 
c/2 \ distance, and hence luminosity scale, of these sources is very interesting. 

This probl em has become pa rticularly important with gamma-ray detection of large numbers of blazars with the 
Fermi LAT (jAbdo et al.l l2009bh . It has been kn own since EGRET tha t radio loud flat spectrum blazars dominate 
I | the b right extragalactic sources in the GeV sky (Hartma nn et al.1 119991 iMattox et al.1 120011 : ISowards-Emmerd et al.1 
^ , 120051) : these objects are fiat-spectrum radio quasars (FSRQ) and radio loud BL Lacs. We have good techniques in 
place to identify likely radio counter parts for the gamm a-ray sources. Extensive spectroscopy campaigns, especially on 
the EGRET-likc 'CGRaBS' sample ([Healev et al.ll2008l ) have provided nearly complete redshifts and characterizations 
£f) . of the FSRQ. However, despite extensive observation, including 8+m telescope integrations, less than half of the BL 
— i ■ Lac counterparts h ave spectroscopic IDs. In some cases, high quality spectroscopy can give useful low er limits on 
the source redshift jSh aw et~alll2009h . but many remain unconstrained. In the first LAT blazar catalog (jAbdo et al.l 
l2009a|) . the problem is even more pronounced since the excellent high energy response of the LAT favors detection of 
hard spectrum sources. BL Lacs are substantially harder ((r) ss 2.0) than FSRQ ((r) sa 2.4) and so provide a larger 
fraction, ~ 40% of the Fermi blazar sample, than was seen by EGRET. 

We attempt here to constrain the redshifts of radio loud BL Lacs which have shown no convincing spectroscopic red- 
shift solution, despite sensitive observations on large telescopes, by high dynamic range imaging searches for the AGN 
host. Our targets are drawn from the 'CRATES' catalog of bright flat-spectrum \b\ > 10° radio sources (jHealev et al.1 
120071 ). Specifically, in this program we targeted the subset of CRATES consisting of i) sources selected as likely 
gamma-ray blazars before the Fermi mission (i.e. CGRaBS sources, Healey et a l. 2008) and ii) sources that are likely 
counterparts to LAT sources detected early in the mission (jAbdo et al.H2009aL f2010h All sources are at Declination 
> —20° and all have been shown to display featureless optical spectra (jHealev et alj 1200a : IShaw et al.l 120091 [2010). 
with sensitive (4m+ class) observations. We were able to image ~ 2/3 of the sources satisfying these critera at the 

time of the observing c ampaigns. 

It has been claimed (jSbarufatti et al.1 120051 ) that BL Lac host galaxies detected with HST imaging are remarkably 
uniform giant ellipticals with Mr — —22.9; accordingly host detections can give redshift estimates and upper limits 
on the host flux can give lower limits on the distance. Here we do not test this assumption or the possible biases 
that would select a modest magnitude range for the detected host sample. We simply apply the method to extract 
redshift constraints, noting that while individual redshifts are doubtless imprecise, the estimates can still be useful for 
statistical studies of BL Lac evolution and as a guide to and comparison with other methods of redshift estimation. 
Conversely, when precision spectroscopic redshifts become available, our measurements can be used to help constrain 
the host evolution. 

Since we will be searching for de Vaucouleurs profile excesses in the wings of the stellar PSF of the bright BL Lac 
core, we require good natural seeing and a moderate field of view for adequate comparison stars. For example, while 
near-IR adaptive optics can deliver superior PSF cores, the wings of the source at >90% encircled energy are extensive 
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and often quite variable over the small corrected FoV. This prevents the accurate PSF modeling and subtraction 
required to obtain host measurements whose integrated magnitude can be 10% or less of the core flux. 

2. OBSERVATIONS & DATA REDUCTION 

All observations were carried out using the WIYN 3.6m telescope at Kitt Peak National Observatory, which features 
good natural seeing and a ~ 10' FoV with the standard mosaic imagers. 

2.1. OPTIC Camera 

We employed OPTIC on the nights of 2008 Oct. 31 - Nov. 2 for our initial imaging run. This 2x2kx4k camera 
covers 9.6' at 0.14"/pixel and features Lincoln Labs orthogonal transfer CCDs, which allow rapid electronic (OT) 
tracking using an appropriate in-field guide star. With a relatively short (25s) read-out time, we were able to combine 
short exposures that avoided core saturation with longer exposures that probe for the host galaxies in the PSF wings. 
The targets were imaged with the SDSS i' filter KPNO 1586. We observed 31 targets, with total on-source time 
ranging from 210s to 2100 s (Table [2]). Sky conditions were photometric through the run, and we were able to employ 
OT guiding for the bulk of the sources, with correction rates up to 50 Hz. The final i' image quality had a median 
FWHM of 0.65"; the best 25% of the sources had image stacks with FWHM < 0.5". A period during which mirror 
cooling was lost resulted in some observations with poor PSF (up to 1"). Software problems in addition cost several 
hours. 

While OT guiding significantly improves the source PSF, which is crucial for our science goals, it does create some 
complications during data reduction. In particular, flat-field frames must be created for each individual exposure, 
with weighted exposure per pixel following the history of OT offsets during each individual exposure. This creates 
some challenges in assembling the final mosaic images, as defects can affect adjacent pixels. The 'smoothing' of the 
flat-field response also affects the noise statistics of the final image. The data were processed using routines in the 
IRAF package mscred. With the excellent PSF (as small as 0.4") achieved for some frames, we required superior 
camera distortion corrections at each position angle; for some orientations a lack of field stars limited the accuracy of 
the astrometric solution. After registration and exposure-weighting, median combined image stacks were prepared for 
further analysis. 

2.2. Mini-Mosaic Camera 

Observations were made using the Mini-Mosaic camera (MiniMo) on the nights of 2009 March 24-25, under highly 
variable conditions. The 2x2kx4k mosaic similarly has a plate scale of 0.14"/pixel and covers a field of 9.6 arcmin. 
We again employed the SDSS i' filter. Because of the long (~3 minute) readout time of the Mini-Mosaic camera, 
most objects had 3-5 dithered exposures of 300s each. These relatively long exposures meant that many of the BL 
Lac targets had saturated cores, despite the typically poorer final image quality during this run (median final FWHM 
0.84"). The exposure sequences and final DIQ for each obje ct are listed in Tab le El In all, 18 BL Lacs were observed. 
Since early LAT detections had been announced by this run (|Abdo et al.ll2009bl) we were able to specially target known 
gamma-ray emitters lacking redshift solutions. The sky transparency was variable and half the run was lost to high 
winds and poor seeing. 

MiniMo employs conventional CCDs, and so standard reductions were conducted with the IRAF package mscred. 
Camera distortions were mapped using USNO A2.0 catalog stars; these were also used to assign a WCS to each frame. 
After processing and bad pixel correction the individual dither frames were aligned using unsaturated stars near the BL 
Lac, exposure weighted and median combined using standard sigma-clipping algorithms to produce our final images. 
MiniMo suffers from 'ghost' images due to amplifier cross-talk. Care was taken during target acquisitions to ensure 
that these did not fall in the vicinity of the BL Lac. We made no attempt to correct these ghosts, but avoid using any 
affected comparison stars. 

3. CALIBRATION 

A number of the target fields were covered by the SDSS, so we were able to establish our zero point directly 
from standard aperture photometry of unsaturated field stars (16 < m*/ < 19). For the OPTIC run only six fields 
had suitable stars. The median rms of these zero point measurements was 0.052 mag. While this was significantly 
larger than the intrinsic SDSS photometric errors, the scatter is smaller than our typical host photometric errors 
and much smaller than the scatter in the claimed absolute magnitude of the BL Lac host (8M = 0.5 mag). For the 
exposures lacking SP SS comparison stars, we applied the mean zero point after correction for atmospheric extinction 
(|Massev et al.l [20021 ) . Given that conditions were photometric throughout, we are conservative in assigning the 0.068 
scatter between the zero point measurements as the final zero-point error. 

While conditions were variable during the MiniMo observations, 12/18 fields had stars suitable for direct SDSS 
cross-calibrations. In particular for four observations logged as cloudy we were able to establish in-field calibration. 
Taking the rms zero point from the remainder of the fields, we established, after extinction correction, a zero point 
with scatter of 0.036 mag. This was applied to the remaining target fields. 

4. PSF MODELING 

Accurate PSF modeling is essential to extract the BL Lac host from the wings of the nuclear point source. We 
built an effective PSF in the stacked image of each BL Lac using the IRAF/DAOPHOT psf routine. PSF stars were 
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Fig. 1. — PSF variation across the field for OPTIC data. Left: FWHM variation with distance from the field center in the image of 
J0110+6805 (observed without OT guiding). Right: variation with offset from the star used for rapid guiding near J0612+4122. 

selected from those that appeared in all sub-frames making up the image stack, had stellar FWHM and appeared 
isolated. PSF stars were selected by hand and were individually checked for consistency with a stellar FWHM. In 
practice the process was iterative; some otherwise desirable PSF template stars had faint unresolved sources in the 
wings, so we modeled and subtracted these sources before recombining for an improved PSF. Point sources as faint 
as i' ~ 25 were removed in this manner. Each template PSF has a radius of 5.6" and covers over 11 FWHM even 
for our worst DIQ images. The background level for each PSF star was determined from the mode of a 2—4" wide 
annulus starting at 7" from the star. In a few of the images with poor seeing the background was determined from 
larger radius, to ensure insignificant contamination from the PSF wings. To best model the PSF wings, we also used 
some stars with nonlinear cores. Above the non-linearity level (30,000 ADU for OPTIC, 40,000 ADU for MiniMo), 
the cores were masked in the PSF stack. 

We found that a PSF model with linear gradients across the image (|Davislll994l ) was required to take full advantage 
of the small FWHM of the best OPTIC data. To guide our fitting we examined the spatial dependence of the FWHM. 
Three terms might be expected to contribute to PSF degradation. First, we expect the efficacy of the low order 
atmospheric corrections to fall off with distance from the guide star. We attempted to mitigate this effect by guiding 
in the target quadrant, whenever possible. Second, distortions in the camera system are expected to increase away 
from the optical axis (assumed to be near the field center) . Finally, since our principal frame registration is on the BL 
Lac core, registration or residual rotation errors could degrade the PSF as a function of distance from the target itself. 
To test for these effects, we collected FWHM measurements of unresolved sources across a number of fields. In no case 
did we find the PSF was primarily correlated with distance to the target either in the individual exposures or in the 
final combined images, implying that stacking error was not important. In some cases the PSF degradation appeared 
correlated with distance from the optical axis, implicating camera distortions, but for other images the distance from 
the guide star appeared to dominate (see Figure [lj . 

With the poorer seeing experienced during our MiniMo run, we found that the spatial variation of the PSF was much 
less important, with the FWHM varying by no more than a few ~ 0.1 pixels (~ 0.03") across the field of view. Again 
stacking errors were negligible. However for consistency we fit for and applied linear gradients in the PSF model. 

Thus, in forming the final model PSF for the analysis, we were able to draw stars from across the fields avoiding 
only the outermost corners of the frames, where there was notable degradation. In general between ~ 1/6 and ~ 1/5 
of the PSF stars used were saturated in their cores (only ~ 1/9 with MiniMo, given the poorer seeing). Tables [2] and 
[3] list the number of PSF stars used in each image. The median number of PSF stars for each field was 30 for both 
OPTIC and MiniMo. 

5. HOST GALAXY MODEL FITTING 

Previous studies (jScarpa et all 120001; lUrrv et alj[2000t lO'Dowd k Urry|[2005l ) indicate that BL Lac hosts are better 
fit by de Vaucouleurs profiles than exponential disks. We therefore adopt simple de Vaucouleurs profiles with fixed 
Sersic index 4, convolving the model with the image PSF, in our galaxy fitting. We did find examples of disturbed 
galaxies and hosts with faint companions, but no host appeared disk-dominated. Our BL Lac image model includes 
up to five parameters, introduced hierarchically. We start by fitt ing for nuclear (PS F) flux plus possible host. We 
assume that the host galaxy isophotes are centered on the nucleus ()Falomo et al.|[2000h . and fit or bound the host flux, 
initially assuming a a fixed effective radius of R c = 1.64" (10 kpc at z = 0.5; we assume ft m = 0.3, Ha = 0.7 and 
Hq = 70 km s _1 Mpc^ 1 throughout this paper). We next allow the de Vaucouleurs profile R c to vary, when required 
by a high significance for a measurement of R e in the fit. Finally for a few BL Lacs, we find a significant detection of 
host ellipticity; in these cases we include the ellipticity e and position angle 9 as additional free parameters in the fit. 
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After fitting the host normalization, we report galaxy magnitudes by extrapolating the integrated light to infinity. 

To prepare for host fitting, we determine the centroid position of the core PSF, using the DAOPHOT task peak, 
which matches to the model PSF, excluding pixels above our (conservative) non-linearity threshold (30,000 ADU for 
OPTIC, 40,000 ADU for MiniMo), and delivering centroid coordinates accurate to ±0.01 pixels. Since a number of our 
images were saturated, this clipping was important for obtaining accurate centroids. In almost all cases, the galaxy 
surface brightness dropped below the background surface brightness fluctuations within our standard 5.6" PSF radius. 
Accordingly the BL Lac image model was fit to 11.2" x 11.2" cutouts from the image stack centered on the BL Lac 
core pixel. DAOPHOT provides a 2x-oversamplcd PSF model. This was re-binned to the original pixel size, using the 
precisely determined centroid and then used in the subsequent fitting. This PSF centering and re-sampling proved to 
be quite accurate - dipole structure was never evident in residuals to our final BL Lac model, whereas with even 0.05 
pixel shifts for the BL Lac core position such residuals were obvious. 

We apply our BL Lac fit to a subset of the pixels in the cutouts around each BL Lac. First all pixels beyond 
5.6" (outside the PSF model) are masked. Pixels associated with resolved objects within the fitting radius but more 
than three FWHM from the BL Lac core are also masked (in practice we exclude rectangular regions, see Figure [2]). 
Unresolved objects more than three FWHM away are fit to the PSF model and then treated as fixed background. 
Occasionally, very bright nearby stellar objects show significant residuals at the PSF core after subtraction; in these 
cases a few pixel rectangular region is masked at the core before the final BL Lac fit. The occasional column or pixel 
defect or cosmic ray are similarly masked. Finally we mask the pixels within 1 FWHM in diameter of the BL Lac core. 
In roughly 25% of the OPTIC images this was required in any case due to non-linear ADU levels. However, we found 
that sampling noise at 1 FWHM could dominate the fit statistics, so a consistent exclusion of pixels at the core was 
used for all objects. For a few of the brightest objects, we extended the core exclusion region to 1.5 FWHM, while for 
a few of the faintest, we used 0.5 FWHM. The final results were not sensitive to the precise value used. 

For each BL Lac model the background was determined from an annulus starting at 7" containing 2-3 x the pixel 
count of the PSF model region; for a few poor seeing images and BL Lacs with very extended hosts we used a larger 
annulus. 




Fig. 2. — Preparing the fitting region for J0814+6431 (image FWHM 0.45"). One resolved companion and pixels near the centroid of 
the nuclear component are masked. The corrected FOV is 7.7" in radius; N up, E left. 



5.1. Noise Model & x 2 Minimization 
For all multi-parameter fits, the parameters are determined by minimization of the quantity: 

r = £^£ m 

i 1 

Where the sum is over the {i} unmasked pixels in the fitting region, Ii is the observed intensity in ADU, M, is the 
model intensity in ADU and of is the expected variance. The model is the sum of the constant background, PSF and 
host galaxy terms, where the analytic host model is integrated over a given pixel to determine the galaxy contribution. 
The variance is estimated according to simple Poisson statistics to be: 



- Put ' "-n (2) 

g x (n im - 1 ) 

Where for OPTIC the gain was g = 1.45 e^ADU^ 1 and for MiniMo 1.4 e _ ADU -1 . The read noise R is R = 4 e~ 
(OPTIC) or 5.5 e~ (MiniMo) and our final image stack is the from the median of ni m normalized single exposures. 



L x q + R 



Imaging Redshifts for BL Lacs 



5 



This is a simplification since for OPTIC the orthogonal transfer during integration correlated the sky noise (likely an 
artifact of imperfect flat-fielding) . Also the residual pixel noise in the unmasked region adds to the variance for sharp 
images. Nevertheless we find that this variance is a reasonable hrst approximation to the local noise, and we can 
determine the best-fit model by % 2 minimization. We use the simulated annealing algorithm from Press et al. (1986). 

The final size of the simplex in N + 1 dimensions gives the relative errors between the N estimates of the fitting 
parameters. We would like to relate this to the la statistical error, to provide a standard statistical error estimate 
for each parameter. We do this by scaling to an externally estimated la error determined for each point source flux 

&PSF- 



We determine this error by analogy with standard aperture photometry. In our case the central FWHM is masked 
in the fit so we compute the fractional error in the PSF amplitude for an annulus outside the masked region of width 



Where N ann is the source flux (ADU) within the photometry annulus, n ann is the number of pixels within the 
annulus, a s k y is the standard deviation of the background (ADU), and n s ky is the number of pixels used to determine 
the background value. The first term under the radical is Poisson noise. The second is a statistical uncertainty due 
to the random fluctuations and unresolved sources in the background. The final term accounts for the possibility of 
a Poisson error in background level. The first term generally dominates, and thus farm is ~ v2 x f ap as the annulus 
contains roughly 1/2 the flux of an aperture without the central mask. We found a typical f a nn for our PSF stars of 
~l-2xlCT 3 . 

With the large number of counts for the point sources measured here, even small centroiding errors can in principle 
contribute to the error, since a varying fraction of the point source falls in the central masked pixels. To test for this 
in each PSF fit, we artificially shifted the model PSF position by ±0.01 pixels, the maximum centroiding error, in 
each coordinate. The rms variation in the included counts was computed to estimate / cen t, the fractional error in flux 
estimation due to imprecise centroiding. Although this error only became comparable to farm for a few of the brightest 
point sources, we sum it in quadrature to compute our statistical photometry error as apsF — Npsf\/ (farm + /cent): 
where Npsf is the total number of counts for the point source in ADU. 

As noted above, the size and shape of the simplex determined from our \ 2 minimization provides relative fitting 
errors. We adjusted our convergence criterion such that the rms spread about the mean fit value for the PSF amplitude 
was always close to apsF which, since the nuclear point source always dominates the model, is a very good estimate of 
the true model-fit statistical error. Thus, by rescaling the simplex amplitude such that the rms spread in the nucleus 
flux estimates was equal to this photometric error apsF, we have normalized the simplex-determined error to be la 
(statistical) in each quantity. We call this error a sta t- 



Since our models are dominated by the bright core PSF, which is different for each exposure, and since we have a 
limited number of stars in each frame to generate a PSF model, our final uncertainty must be dominated by sys tematic 
errors due to imperfections in this model. This is in contrast to the HST- based fitting of iScarpa et al.l (|2000| ). where 
the PSF is, of course, stable and well-modeled and unresolved background structure dominated the final uncertainty. 
We therefore wish to estimate our systematic PSF uncertainty. 

We do so by fitting for a de Vaucouleurs host around unresolved stellar sources. In each image stack we choose such 
test stars relatively close to the BL Lac with comparable brightness. On average we were able to measure seven such 
stars per image. Ideally none would have been used in forming the PSF, but since these are the stars most suitable 
for PSF modeling, they were largely included. Before fitting, each star was treated exactly as its associated BL Lac, 
including the same masking of central pixels. We then fit for the PSF amplitude and 'host' amplitude, with the de 
Vaucouleurs R c = 1.64" fixed at our default radius. The fit host amplitude could either be negative (indicating excess 
wings on the model PSF) or positive (PSF model too narrow). Note that decreasing the host angular size would 
weaken constraints on its amplitude, since it would be more highly covariant with the PSF. However, our adopted 
value (appropriate for a typical R e = 10 kpc at z = 0.5) is reasonably conservative; even at the minimum angular size 
at z ~ 1.6, this would only decrease by 28%. 

The dispersion in the amplitudes of the fitted 'host' fluxes provides an estimate of our systematic uncertainty. 
However, since the stellar fluxes covered a substantial range, we elected to scale the 'host' flux significance to the 
statistical error on this flux for each test star <7 s t a t The rms spread in the 'host' fluxes fit around the test stars, in 
units of Cstat, then provides an estimate of the systematic error, for each individual BL Lac field. Figure [3] shows 
the distribution of 'host' significance (iVhost/Vstat) for the test stars in the OPTIC and MiniMo fields, along with the 
distribution of the rms significance for the individual fields. 

The rms significance of these 'host' fits for OPTIC and MiniMo were 8.3a s t a t and 8.0a s t a t, respectively. This suggests 
that our statistical (Poisson) errors on the host flux substantially underestimate the true errors by ~ 8x. However, it 
should be noted that the distributions are approximately normal and that the mean fitted host flux is not significant. 
For the OPTIC data, we obtain — 0.7 a s tat while for MiniMo the average is 0.1a s t a t- In view of the systematic errors we 
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Fig. 3. — Statistical significance of 'host' detections around field stars. Left: OPTIC stars for 31 fields, Right: MiniMo stars for 18 fields. 
The dark histograms show the distribution of the rms significances for the individual BL Lac fields. There is a modest tail, but most cluster 
near the global rms. These individual rms values give our systematic error estimates for the individual BL Lac fields. 

expect random offsets of ±0.6cr stat and ±0.7a sta t for the two sets of test stars. We conclude that the frame-to-frame 
spread in upper limits on surrounding host galaxies are appreciably larger than expected from pure Poisson statistics 
- this can be attributed to errors in forming the individual PSF models and to unresolved structure in the PSF stars' 
and test stars' backgrounds. However there is no overall introduced host flux; implying that our PSF model is a fair, 
albeit uncertain, representation of the data. 

To be conservative, we adopt for our final error the sum of the statistical and systematic errors, namely <7fi n = 
c S iat(l + °~syst)i where a sys t = (Nhost/ o~stat)rms is the rms value of the individual test stars in the stacked image of 
an individual BL Lac. This is conservative since while we expect the statistical errors to have a normal distribution, 
it would be very surprising if the much larger systematic errors had wings as large as a Gaussian out to many times 
the FWHM of the distribution. We tested this error estimate by inserting Poisson realizations of model hosts around 
isolated stars in several fields. For artificial galaxies assigned the total counts at the 3tJ/i„ detection limit, we recovered 
host fluxes with rms errors of 0.7tx/i„. In no case was the difference between the injected and recovered counts larger 
than lcr fi n . 

5.4. Final BL Lac Fitting Procedure 

We constrain the properties of the BL Lac host by the hierarchical fit. We first fit for a circularly symmetric host 
with fixed angular radius R c = 1.64", just as for the test star measurements. The final uncertainty on the host counts 
is taken to be afi n where the statistical error is taken from the simplex rescaled to the aperture photometry error for 
the BL Lac core, while the systematic error is the multiple of this statistical error determined by 'host' fits to the 
test stars in this BL Lac's field, as described above. If the host flux is larger than 3cr/i„ (which is dominated by the 
systematic uncertainty), we deem it significant. If the significance does not reach 3<7/i n , we infer an upper limit of 
Nhost + 3cr f in , or 3cr/ irl , whichever is larger. 

If there is a significant detection of a BL Lac host, we next re-run the fits allowing the host angular radius R e to 
vary. The only fitting constraint is R c > 0. If the best fit value of R e is less than the FWHM for the image, we do 
not deem it significant. A value for R c is only quoted if R e > 3<r/i„, where the final statistical plus systematic error is 
estimated from the re-scaled simplex, as above. Significant estimates of R e were found for 25% of the OPTIC-measured 
BL Lacs and 33% of those measured with MiniMo. Note that the overall flux and its fit error were free to vary, as 
well. As a consistency check, we note that in no case where a significant R c measurement was found, did the host flux 
measurement decrease in significance to below 3cr. 

Finally, we attempted to fit for a significant host ellipticity, with two additional parameters e and the position angle 
9 added to the model. These parameters both exceeded 3cr/i„ significance for only two of the BL Lacs observed with 
OPTIC. Again the host flux and size R e were free in this final fit. 

Tables [5] and |31 contain the final fit nucleus and host amplitudes (converted to magnitudes, see below) and the host 
physical size (when available). Details of the observations (exposure lengths, FWHM of the stellar sources and the 
number of stars used in the PSF model) are also listed. 

6. REDSHIFT ESTIMATES & LOWER LIMITS 

Since iSbarufatti et al.l (|2005| ) have argued that BL Lac host galaxies are standard candles with Mr = —22.9 ± 0.5, 
we can use our measurements and upper limits on the host flux to extract redshift estimates and lower bounds. To do 
th is we compute an i' Hu bble diagram for our KPNO 1586 filter. We improve on the similar R-band Hubble diagram 
of ISbarufatti et ail ()2005l ) . by including host evolution, assuming a host (elliptical) formation redshift of Zf orm — 2 
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Fig. 4. — Filter transmission curves used in this study. Solid i' , 
dashed R, dotted HST F702W. We convert from the assumed host Fig. 5. — Hubble diagrams calculated in KPNO w004 R and 
z = R magnitude to our observed i' magnitude. Note that the i' K PNO 1586 i' . The an alytic estimate for the R Hubble diagram 
cut-on is redder than used in previous studies. of [Sbarufatti ct al. (2005) is shown for comparison. 

(jO'Dowd fc Urrvl l2005). With this assumption, we adopt an elli ptical galaxy spectrum computed at the appropriate 
age from the PEGASE model of lFioc fc Rocca-Volmerangel ()1997f ). whe re the age is computed for each 2 obs < zj orm for 
our standard cosmology. To lock the overall normalization to the Mr of lSbarufatti et alT(|2005[ ). we fixed the flux of the 
z = 0, evolved elliptical spectrum to Mr = —22.9. This required folding through an R filter. Unfortunately the precise 
KPNO K ron- Cousins R filter u sed is not recorded, except in a private 1995 communication between J. Holtzman and 
Landolt (Holt zman et al.11199 5). Accordingly we use the transmission curve for Kron-Cousins R filter KPNO w004, 
first obtained in 1994; this may be the precise filter used in the original study. In any case, after normalizing with 
this filter at z = 0, we take the model flux from each observation redshift, convert to the observed frame using the 
luminosity distance dh from our standard cosmology and convolve with the i' KPNO 1586 transmission (Figure 2]) to 
obtain the apparent magnitude. 

The resulting Hubble diagram is shown in Figure [5] For co mparison we hav e also computed the i?-band Hubble 
diagram for this evolving model. The curve is quite close to Sb arufatti et al.l (|2005|) 's analytic approximation, al- 
though departures caused by the evolution are evident. Hubble diagram computations with other evolving galaxy 
models yielded similar curves, with differences small compared to those introduced by uncertainty in the host absolute 
magnitude. 

To compar e with our measurem ents, we conv erted our model fit h ost counts (in ADU) to i' magnitudes, including 
atmospheric (Massev et al.l l2002 s ) and Galactic (Schle gei et al.|[T998f) extinction corrections. The establishment of the 
flux scale is discussed above (§3), while the assumed Galactic extinctions from the dust maps at the BL Lac locations 
are listed in Tables [2] and [3] 

Finally, these host magnitudes (and lower limits) are converted to redshifts using our Hubble diagram. Of course 
the host absolute magnitude uncertainty translates into a range of possible redshifts. There are two contributions 
to the full redshift uncertainty. The first is host measurement uncertainty, including both statistical and systematic 
photometry errors. We use the ±cr/i„ range for the maximum and minimum host magnitude to infer asymmetric 
error bars on the host redshifts. This is augmented by adding, in quadrature, the photometry errors estimated for our 
zero point - this increase is however very small. These values are listed as the first error flags in the final columns of 
Tables [2] and In addition, the claimed SR = 0.5 mag dispersion in the host luminosity leads to a redshift range. We 
propagate the z — R absolute magnitude uncertainty through observed i' and hence to Az as a function of redshift 
(see Figure [S]). This nearly always completely dominates the redshift uncertainty and is listed as the second error flag 
in our tables. 

For objects with a fit angular R c , the central redshift estimate is used to convert to proper kpc; the values are 
comfortably close to the otherwise assumed R c = 10 kpc. Finally, for undetected hosts, we take the minimum z allowed 
by the 3<r/i ra host magnitude limit, the photometric scale errors, and the standard value for the host luminosity to 
obtain a lower limit on the host redshift. Again this is listed in the last column of Tables [2] and 

We have treated our measurement errors conservatively, having summed the systematic and statistical uncertainties 
in our error flags and requiring 3 x a fi n detections before claiming a measurement is significant. On the other hand, we 
consider only the assumed dispersion in the host luminosity. Clearly an improved calibration of the host magnitudes 
(if uniform) can make these measurements much more useful. 



7. RESULTS AND DISCUSSION 
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Fig. 6. — Upper and lower error error bars as a function of redshift, inferred from the ±0.5 dispersion in host absolute magnitude Mji at 

2 = 0. 

We were able to detect hosts for 14/31 OPTIC-observed BL Lacs and 11 of 18 observed with MiniMo. The redshift 
estimates varied from ~ 0.20—0.95 and the lower limits from ~ 0.25—0.9. For OPTI C the median value w as z — 0.54 for 
MiniMo z — 0.53. It is interesting to compare this with the HST host detections of IScarpa et al.l ([2000); in this study 
nearly all BL Lacs at z < 0.5 were resolved but only 6/23 sources with known z > 0.5 yielded host detections. While at 
first sight it might seem surprising that a ground-based program could be competitive, it should be remembered that 
the HST snapshot BL Lac survey observed with the F606W and F702W filters. The i' filter blue cut-on is ~ 100 nm 
redder than for the bulk of the HST observations, and so we do not expect to suffer appreciable surface brightness 
attenuation from the 4000A break until z ~ 0.75. Further the typical HST snapshot exposures were ~ 500 s with a 
2.4m aperture vs. 900— 2100s with a 3.6 m aperture. The combination of redder band and deeper exposure allow some 
additional sensitivity to the low surface brightness host wings. Of course, our sources were selected from those lacking 
spectroscopic redshifts, making it likely that these are a higher redshift sub-sample. However it is also worth noting 
that the typical ~ 10 — 30 x core host ratio plotted in Figure [7] represents an even stronger core dominance, compared 
to the HST sample, than one might think since the continua of BL Lacs is quite blue compared to the hosts, and our 
i' bandpass is redder than that of previous R or FQ0QW/F702W studies. 

For 14 of the 25 detected hosts we obtained estimates of the host angular size. The inferred median effective radii 
were 9.2 kpc (OPTIC) and 8.0 kpc (MiniM o) in good agreement with previous BL Lac host estimates (jScarpa et al.l 
120001 ; ItJrrv et al.ll2000l ; iFalomo et aLll2000[) . The other detected hosts had R c measurements with < 3aj in significance 
or best fit radii smaller than the FWHM and thus were judged to be insignificantly resolved. 

In two cases we have a formal detection of ellipticity. For J0211+1051 the small value e = 0.21 accords well with 
the visual impression (Figure ITT]) . In the other case, J0348— 1610, the ellipticity is also small e = 0.14, but here the 
result is suspect, as there is a bright companion within 2R C . Indeed it may be no coincidence that the two objects with 
significant e both have a large number of companions. While the upper limits on e are in some cases constraining, we 
conclude that with our modest ground-based resolution we cannot in most cases probe the host structure. 

One arena where the comparison with HST-based measurements should be fairly robust is the host -nucleus flux 
ratio. Here we find a marked difference with the HST sample. Of the 69 BL Lac hosts resolved by IScarpa et al.l 
(2000), 54% of these had f host / f nucleus > 1. Even counting the 42 unresolved BL Lacs 34% had a nuclear- host ratio 
> 1. In contrast, we find only 1 object of our 49 BL Lacs has f host/ f nucleus > 1- Indeed 63% of our hosts are > lOx 
fainter than the nucleus (Figure [7]). Two selection effects may account for this difference. First, these objects have 
been selected as bright flat spectrum radio core sources. Further a significant fraction of the MiniMo targets were 
in addition known to be active gamma-ray emitters. This implies that the Earth line-of-sight is even more nuclear 
component dominated (on-axis) than for the typical BL Lac. Second, these are the subset of these systems lacking 
previous spectroscopic redshifts. Again, this suggests that they should be even more core continuum dominated than 
the typical BL Lac. 



7.1. Radial Profile Plots 
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Fig. 7. — Host magnitude vs. nucleus magnitude. The dashed lines show hosts lOx and lOOx fainter than the BL Lac nucleus. 
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Fig. 8.— OPTIC radial profiles for fit BL Lacs: Left: J07I2+5033, resolved at z ^ 0.47. Right: J2305+8242, resolved at z si 0.62. 
The dashed lines show the PSF, the dot-dashed lines the host and the solid line the total model. Error flags on individual radial bins are 
statistical only. 

We present here a sample of azimuthally averaged i' instrumental surface brightness plots, compared with our best 
fit model components (Figures 151 — [TU)) . The surface brightness data are measured from excess counts above the above 
the fitted background level measured in annuli about the core position. The model curves are from integrations over 
these radial bins. The error bars shown are la, calculated in a fashion consistent with the variance in §5. We show 
two BL Lacs with host detections using the OPTIC camera, as well as two fits providing upper limits. For MiniMo 
we show one detected host and one unresolved source. 

7.2. Near Environments 

It has been noted (|Scarpa et al.ll2000t lUrrv et aLll2000t iFalomo et al. I2000f) that BL Lac hosts are surrounded by a 
significant excess of nearby galactic objects. Following lO'Dowd fc Urrvl "( 2005( ) . we count companions within a projected 
proper distance of 50 kpc of the BL Lac centroid. When we do not have a z estimate, we count only companions within 
5.9" of the nucleus (50 kpc at the minimum angular diameter distance z ~ 1.6). We count only resolved objects. 
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Fig. 9.— OPTIC radial profiles for fit BL Lacs with host non-detections: Left: J0049+0237, unresolved z > 0.84. Right: J0607+4739, 
unresolved z > 0.54. 
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Fig. 10.— MiniMo radial profiles for fit BL Lacs. Left J1008+0621, resolved at z fa 0.59. Right: J0953-0840, unresolved 2 > 0.40. 



We find at l east o ne compani on for 51% of ou r BL L acs (25/49), in reasonable agreement with the detection rate 
of lUrrv et ail (pOOTl (47%) and lO'Dowd fc Urrvl (pOOl (62.5%). Of course, with our ground-based imaging, some of 
the unresolved nearby sources may represent additional compact galactic companions. The number of companions N c 
within a projected distance of rj_ = 50 kpc, the magnitudes of the two brightest companions (and the comparison of the 
brightest to the host flux Ami, host) and a flag indicating evidence for interaction are given in Table[T] As noted above, 
our two hosts with nominal detections of ellipticity have both large numbers of nearby companions and morphological 
evidence for interaction. These may in fact represent recent mergers in compact groups. J1440+0610 has an extended 
companion with centroid at r± > 50 kpc which extends within 50 kpc. For two objects (J0050— 0929 and J0743+1714) 
the companion at small radius may be physically contiguous with the undetected host. Clear evidence for interaction 
was seen only for the better-seeing OPTIC data, so additional companion interaction is likely. 



7.3. Comparison with Other Redshift Estimates 

As a part of our on-going spectroscopic campaign to obtain identifications and redshifts for the -Ferrm-detected 
blazars, some of these imaging targets have had additional spectroscopic exposure. In a few cases we obtain true 
spectroscopic redshifts z s . In several other cases we have lower limits on the redshift z a bs from detection of intervening 
intergalactic absorption-line systems. Finally, in several cases we were able to use spectroscopic limits on absorption 
line-strengths (typically Ca H & K and G-band limits), together with the assumpt ion of a uniform host magnitude (as 
in this paper) to place lower redshift limits zhk on the BL Lac (jShaw eEaDHQOl). 

We have obtain ed a direct redshift measurement for only one of these targets, using Keck LRIS spectroscopy 
(Sh aw et al.ll2010[ ). The value is within lcr (statistical+systematic) of the imaging estimate: 

J2022+7611 zr = 0.49±g:g|±g:gl (z s = 0.584). 
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TABLE 1 

List of brightest observed companions within r± = 50kpc. Except for those marked f, the brightest companion is also the 
closest. We classify the morphological evidence for interaction as obvious "Y," possible "M" or absent "N." 




Fig. 11. — Two BL Lac fields after core PSF subtraction. Left: J0211+1051 at z RJ 0.20, showing the highly resolved host and several 
fainter companions. The host has significant cllipticity; the fit _R C ellipse is shown. Right: J2241+4120 at z fa 0.52, showing the faint wings 
of the host and two fainter companions. No significant ellipticity is detected. 

For t hree objects our imaging zi estimates are in agreement with lower limits from our own spectroscopy (|Shaw et al.l 
2009). The low redshift of J1427+2347 is particularly interesting in view of the recent VERITAS TeV detection of 
this source (Ong et al. 2009). 

J0712+5033 zi = OA7toM2-os>l ( z hk > 0.47). 
J1427+2347 zj = 0.23+°;°^° ^(z ff/f > 0.03). 
J1440+0610 zi = 0.55±8;g±g;Jg(z o6 , > 0.316). 

For two objects our imaging redshift constraints are stronger than those obtained from the spectroscopic constraint 
on the HK line strengths: 

J0909+0200 zi > 0.83 (z HK > 0.54 IShaw et al.l (t2009h. 
J0049+0237 zi > 0.84 (z HK > 0.82 lSbarufatti et al.l (|2006l ). 

For t wo objects our lower limit on the redshift from imaging is not as strong as that from our spectroscopy (|Shaw et al.l 
120091) : 

J0050-0929 zi > 0.27 (z HK > 0.44) and 
J2050+0407 zi > 0.74 (z abs > 0.819). 
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Finally for two objects the imaging redshifts are in disagreement with the lower limits from limits on the HK absorption 
line strength: 

J1253+5301 zj = 0.59±g;g|±g:io ( z hk > 0.77; lAa higher) and 
J1542+6129 zj = 0.39±o;g|±o.°f (z HK > 0.63; 2.2ct higher). 

These last two disagreements may partly be attributed to differences in the host size and slit losses assumed in the 
spectroscopic method. However, statistically we should expect some Zi estimates to be low, as Malmquist bias will 
assure that fluctuations causing detections make the sources appear artificially close. Of course it is also likely that 
outliers in the BL Lac host magnitude distribution exist: these will cause larger disagreements. As we continue to 
collect spectroscopic redshifts and limits for these sources, w e should be able to pro be the fraction of sub-luminous 
hosts and test the utility of the standard candle hypothesis of iSbarufatti et all (|2005l ). 

8. CONCLUSIONS 

We have detected hosts for half of the BL Lac objects imaged in this WIYN campaign. These detections, assuming 
a standard host luminosity, give redshift estimates with a median value of z me d = 0.51. The upper limits on host flux 
for the remaining objects also give reasonable redshift constraints with a median value for the bound of z mo d > 0.61. 
This may be compared with the spectroscopic redshifts already obtained for the two parent BL Lac populations: for 
the BL Lacs in the early Fermi blazar list one finds z mo d = 0.33, while for the radio selected CGRaBS BL Lacs, we 
have z mc d = 0.45 (Figure [T2")) . We can further quantify the difference by making Kolmogorov- Smirnoff comparisons of 
the distributions. While the LBAS and CGRaBS sets are quite consistent, with a KS probability of 0.58, the imaging 
set differs from both with a KS probability of similarity of only 1 — 5 x 10~ 3 . If we include the lower z bounds in the 
distribution the probability drops to < 2 x 1CP 5 . 
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Fig. 12. — BL Lac redshift distributions. Estimates from host detections (filled histogram), and with host lower limits (line histogram). 
For comparison we also show the spectroscopic BL Lac redsh ift distribution from the two parent populations providing our target lists: 
Squares - the CG RaBS fiat spectrum radio loud BL Lacs (Healey et al. 2008) and Circles — the early Fermi gamma-ray-detected BL Lacs 
flAbdo et al.ll2009al) . 

The evident lack of low redshift imaging detections is doubtless a selection effect: such objects will in general show 
the strongest absorption features from the host and thus are most likely to provide spectroscopic redshifts. With 
the ground-based imaging results, we obtain a number of additional redshift estimates, with good efficiency for host 
detection out to z ~ 0.65. The distribution of our upper limits suggests that space-based imaging searches, using near- 
IR filters, can be productive to even higher redshift. In any event, it appears that it would be a mistake to assume that 
the objects in the radio and gamma-ray samples lacking spectroscopic redshifts are similar to those with spectroscopic 
distances. The imaging results suggest, in contrast, that these other sources are appreciably more distant and, on 
average have higher luminosity. It will be important to include these imaging z estimates and bounds in attempts to 
measure the BL Lac luminosity function. 

One other major conclusion of this analysis is that host searches should be at the i band or redder. Here the 
contrast with the relatively blue nuclear continuum flux is improved: we have seen that this allows host detections in 
sources with extreme core dominance. The increased sensitivity to the low surface brightness halos of the redshifted 
hosts should also materially improve the detection fraction at z > 0.5. Again this will be important for study of the 
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population and evolution of BL Lac hosts galaxies. Finally, it will be important to pursue spectroscopic confirmations 
of these redshift estimates, as only with such data will it be possible to test and extend the claimed uniformity of the 
BL Lac host population and to probe its possible evolution. 
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We thank the staff of the WIYN observatory for excellent support and in particular Steve Howell for help with OPTIC 
and Di Harmer for the MiniMo set-up. We in addition thank J. Windschitl for sharing notes on OPTIC reduction 
techniques and Mike Shaw and Steve Healey for help with comparison to the spectroscopic BL Lac sample. 

APPENDIX 

Here we note peculiarities of the fields of individual BL Lacs, particularly those that required amendment to the 
standard analysis described above. 

J0050 0929: This high Galactic latitude field is particularly sparse, with a deficit of suitable of PSF stars (13). The 
BL Lac host shows asymmetric structure to E, which we flag in Table []] as due to a close companion. Even with a 
relatively good 0.59" FWHM its nature is not clear. In fact, the fixed R c fit suggests a resolved host with significance 
~ 2.5cr/m, but gives an exceptionally high x 2 /D.O.F. ~ 12. This implies that the source is not well modeled by a 
de Vauco uleurs host. Thus we interpret the 3<r non-detection as a redshift lower limit, inferring a distance consistent 
with the dShaw et al.ll2009l ) spectroscopically derived bound z > 0.44. 

J0110+6805: This target is at |6| < 10° and so is not a member of the CRATES catalog. It is otherwise known as 
4C+67.04. An unresolved object of comparable brightness appears ~ 2 FWHM from the AGN nucleus. We did a PSF 
subtraction of this star and masked the peak pixels. Since our PSF fits are sequential, not simultaneous, the substantial 
overlap with the BL Lac core PSF introduced additional uncertainty in the BL Lac host fit. We approximate this by 
adding the aperture photometry PSF uncertainties of the neighbor star and BL Lac nucleus in quadrature, using this 
as the nucleus statistical uncertainty and propagating through the rest of the analysis. 

J0203+7232: A bright i — 14 star appears ~ 4" from the BL Lac centroid. As for J0110+6805 we increase the 
nucleus photometric error by adding in the uncertainty in the amplitude fit for this star. 

J0211+1051: The detectable BL Lac host extends beyond the standard fitting radius 5.6"; accordingly we extend 
the radius of the fitting region to 11.2", and increase the size of the sky annulus. 

J0348 1610: A very bright i ~ 11 field star 20" from the BL Lac produced noticeable scattering wings across the 
fitting region. We estimated this scattered flux by radially averaging the stellar wings over regions free of background 
objects and then subtracted this flux from the fitting region to remove the small induced gradient. After subtraction, 
we applied the standard BL Lac fitting. 

J0743+1714: Two bright (i ~ 14.5) stars are ~ 20" and 25" from the BL Lac. Subtraction of a radially averaged 
template of the scattering wings, as for J0348— 1610, produced a flat background acceptable for the BL Lac fitting. 

J0814+6431: This BL Lac extends beyond the standard fitting region. This is increased to 7.7"; the sky annulus 
radius is increased accordingly. 

J0817 0933: Two bright (i ~ 14) stars are ~ 25" and 30" from the BL Lac. Subtraction of a radially averaged 
template of the scattering wings, as for J0348— 1610, produced a flat background acceptable for the BL Lac fitting. 

J0907 2026: With a bright nucleus and relatively poor seeing, we started the background sky annulus at 10.5". The 
BL Lac fitting still employed the standard 5.6" radius. 

J0915+2933: The scattering wing of a bright i ~ 13 star ~ 50" away was removed, as for J0348— 1610. 

J1253+5301: The scattering wing of a bright i ~ 13.5 star ~ 40" away was removed, as for J0348— 1610. 

J1427+2347: The bright BL Lac nucleus is highly saturated, resulting in bleeding. In this case the central exclusion 
was increased to diameter of 3 FWHM and the bleed trail was carefully masked. The test stars were fit with a matching 
central exclusion region. While the sky background annulus radius was increased to 16.8", the BL Lac model fit was 
performed in the standard 5.6" radius region. 
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u - zu -.01-.04 


J0219-1842 


[180,180,180,180,180] 


0.76 


-68.1 


0.07 


20 


17.65 


19.85^;^ 


0.14 


- 


„ ~~4- 04+ 11 
60 

u - DU -.03-.10 


J0348-1610 


[180,180,180,180,180] 


0.70 


-47.6 


0.09 


24 


16.67 


18.651-Q7 


0.16 


11.1 ± 2.33 


„ L 01 -U OS 

u,o3 -.01-.07 


J0433+2905 


[180,180,300,300,300] 


0.46 


-12.6 


1.49 


30 


17.53 


> 19.20 


< 0.21 


- 


> 0.48 


J0502+1338 


[300,300,300,300,300] 


0.49 


-16.8 


0.99 


28 


18.18 


I9.05i;°» 


0.44 


8.87 ± 1.73 


n 4 t:+.01+.09 
u '^°-.01-.08 


J0509+0541 


[60,60,60,60,60] 


0.53 


-19.6 


0.21 


28 


14.75 


> 18.60 


< 0.03 


- 


> 0.38 


J0527+0331 


[300,300,300,300,300,300,300] 


0.67 


-16.9 


0.31 


27 


18.64 


21.241^ 


0.10 


4.91 ± 1.63 


n Qn+ 05+.14 
u ' au -.04-.ll 


J0607+4739 


[180,180,180,180,180,180] 


0.78 


12.9 


0.36 


30 


15.53 


> 19.56 


< 0.02 


- 


> 0.54 


J0610-1847 


[180,180,180] 


0.83 


-17.3 


0.17 


24 


17.45 


2O.O7I31 


0.09 


- 


n 64 + - 10+ - 12 

U-D4_.06-.10 


J0612+4122 


[130,90,120,120,120] 


0.50 


10.9 


0.39 


28 


16.67 


> 20.28 


< 0.04 


- 


> 0.69 


J0625+4440 


[180,180,180,180,180] 


0.60 


14.4 


0.29 


23 


17.12 


> 20.60 


< 0.04 


- 


> 0.77 


J0712+5033 


[240,120,120,120] 


0.60 


23.9 


0.13 


28 


16.42 


i9.n+-\l 


0.08 


8.68 ± 1.83 


n 47+ 03+. 09 

-.02-.08 


J0814+6431 


[90,90,90,90,90] 


0.45 


33.2 


0.12 


26 


15.51 


18.38±;£g 


0.07 


11.4 ±3.06 


n ocr+.oi+.OT 

u " : ' o -.01-.06 


J0909+0200 


[180,180,180,180,180] 


0.52 


31.4 


0.06 


28 


17.90 


> 20.91 


< 0.06 


- 


> 0.83 


J1813+0615 


[180,180,180] 


0.83 


11.3 


0.41 


106 


17.40 


> 20.31 


< 0.07 




> 0.70 


J190311+554044 


[120,120,120,120] 


0.53 


20.5 


0.12 


36 


16.27 


> 19.75 


< 0.04 




> 0.58 


J1927+6117 


[120,120,120,120,120] 


0.69 


19.5 


0.13 


70 


16.95 


19.57l;J 


0.09 




54 +. 02+. 10 

u, ° -.02-.09 


J2009+7229 


[180,180,180,180,180] 


0.66 


20.2 


0.53 


30 


17.23 


> 20.50 


< 0.05 




> 0.74 


J2022+7611 


[180,180] 


0.64 


21.1 


0.47 


37 


16.68 


19.26lig 


0.09 




4O+ 04+.09 
u, * a -.03-.08 


J2050+0407 


[180,180,180,180,180] 


0.83 


-24.0 


0.18 


29 


17.93 


> 20.48 


< 0.10 




> 0.74 


J2200+2137 


[300,300,300,300,300] 


0.67 


-26.0 


0.20 


32 


18.50 


> 20.64 


< 0.14 




> 0.77 


J2241+4120 


[300,300,300,300,300] 


0.70 


-15.2 


0.47 


31 


17.68 


19.46± 11 


0.19 


9.60 ± 1.00 


n .,,+.02+. 10 
U ' OZ -.02-.09 


J224356+202101 


[90,60,60] 


0.56 


-33.4 


0.09 


39 


15.10 


> 18.64 


< 0.04 




> 0.39 


J2305+8242 


[300,300,300] 


0.94 


20.6 


0.46 


31 


19.99 


19.971 '\b 


1.03 




62 + 04+ - 12 

u,u -.03-.10 


J2346+8007 


[300,300,300,300] 


1.01 


17.6 


0.45 


68 


17.90 


> 21.07 


< 0.14 




> 0.87 



a J0110±6805=4C±67.04 

TABL E 2 

Summary of OPTIC results. 4+4 names are from CGRaBS I|Healey et AL.II2008T 1. 6+6 names are from CRATES (|Healey et al.|[2007Ti . There are also two > 3<r 

DECTECTIONS OF HOST ELLIPTICITY: J0211 + 1051 WITH e = 0.21 ± 0.02 AT 6 = 19° ± 10° AND J0348-1610 WITH e = 0.14 ± 0.04 AT = 32° ± 15° WITH ANGLES MEASURED N 

THROUGH E. 



Name 


Exposures (s) 


Seeing (") 


6(°) 


A 7 


NpSF stars 


Nucleus 


^'host 


J host 1 /nucleus 


Re (kpc) 


z 




J024330+712012 


[300,300] 


1.09 


10.4 


1.50 


51 


16.95 


> 18.70 


< 0.20 


- 


> 0.40 




J065047+250304 


[60,300,300,300,300] 


0.68 


11.0 


0.19 


31 


16.14 


19.031 h, 


0.07 




0.45±;°2t 


.09 

.08 


J0743+1714 


[300,300,300] 


0.65 


19.0 


0.07 


11 


18.43 


> 20.82 


< 0.11 




> 0.81 


J0817-0933 


[300,300,300,300] 


0.78 


14.3 


0.20 


42 


16.61 


20.371'ig 


0.03 


12.5 ±0.58 


n 71 +.05+ 

u - (i -.04- 


.11 
.11 


J083543+093659 


[300,300,300] 


0.83 


27.5 


0.10 


30 


17.93 


> 18.89 


< 0.41 




> 0.42 




J0907-2026 


[300,300,300] 


0.97 


18.0 


0.34 


35 


14.46 


18.32+'2o 


0.03 


8.05 ±0.39 


0.34l;°il 


.07 
.0(1 


J091553+293326 


[300,300,300] 


0.68 


42.9 


0.05 


23 


15.71 


18.451^ 


0.08 


8.50 ±0.16 


0.36±;°|1 


.07 
.06 


J095301-084034 


[300,300,300] 


0.78 


33.9 


0.92 


34 


15.29 


> 18.71 


< 0.04 




> 0.40 


J1008+0621 


[300,300,300] 


0.73 


46.0 


0.05 


26 


17.68 


19.801 a 


0.14 


14.8 ±4.32 


0.59±-.g|± 


.11 
.10 


J103742+571158 


[300,300,300,300,300] 


0.76 


51.8 


0.01 


20 


15.97 


> 19.95 


< 0.03 




> 0.62 




J105912-113424 


[300,300,300] 


0.85 


42.7 


0.05 


21 


15.86 


> 19.08 


< 0.05 




> 0.46 




J124313+362755 


[300,300,300,300] 


0.86 


80.5 


0.02 


26 


15.72 


19.32±;SS 


0.04 




0.50l-g*l 


.10 

.on 


J125311+530113 


[300,300,300,300,300] 


0.93 


64.1 


0.02 


22 


16.83 


19.831;°* 


0.06 


7.37 ±0.73 


0.59 +: ^+ 


.11 

.10 


J142700+234802 


[300,300,300] 


0.71 


68.2 


0.11 


28 


13.99 


17.301;°° 


0.05 




0.23+'^+ 


.or, 

.04 


J144052+061033 


[300,300] 


0.98 


56.6 


0.07 


26 


16.89 


19 - 59 ±:io 


0.08 


7.24 ±2.31 


0.55+'f 2 + 


.10 
.09 


J154225+612950 


[300,300,300] 


1.02 


45.4 


0.03 


33 


14.98 


18.681;^ 


0.03 




0.39+'^+ 


.OS 
.07 


J1624+5652 


[300,300,300] 


0.89 


42.3 


0.02 


33 


17.93 


20.07l'i° 


0.14 




0.64± : °f± 


.12 

.10 


J1749+4321 


[300,300,300] 


1.09 


29.2 


0.07 


30 


17.41 


> 20.73 


< 0.05 




> 0.79 



TABL E 3 

Summary of MiniMo results. 4+4 names are prom CGRaBS i|Healey et al.|[2008T ). 6±6 names are from CRATES (|Healey et al.|[2Q07Ti . No {e, 0} was found with 

SIGNIFICANCE > 3tT FOR ANY BL LAC IN THESE DATA. 



